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"Electromagnetic  Interaction,  Thermal  and  Mass  Transfer  Modeling  of  the 
Photothermal  Modulation  of  Mie  Scattering  (PMMS)  Spectroscopy  of  Aerosols" 

John  P.  Barton 

University  of  Nebraska-Lincoln 
January  15, 2000 


I.  Introduction 

The  procedure  of  photothermal  modulation  of  Mie  scattering  (PMMS)  spectro^opy  is  a 
promising  technique  that  can  potentially  provide  the  sensitive  in  situ  measurement  of  unportant 
aerosol  characterization  parameters  such  as  particle  size  distribution,  index-of-refraction,  and 
composition.  In  the  PMMS  procedure,  the  intensity  of  an  infrared  laser  beam  (e.g.,  a  CO2  laser) 
is  modulated  in  time  (i.e.,  chopped)  and  passed  through  the  aerosol  cloud.  The  dianieters  of  the 
droplets  in  the  aerosol  cloud  subsequently  also  modulate  in  time  due  the  associated  time-vai^g 
heating,  evaporation,  and  recondensation.  A  second  laser  beam,  in  this  case  a  visible  light  laser 
beam  (e.g.,  a  He-Ne  laser),  is  passed  through  the  aerosol  cloud  and  the  effect  of  the  tme- 
variation  of  the  droplet  diameters  due  to  the  modulated  infrared  laser  heatmg  is  observed  by 
monitoring  the  time-variation  of  the  elastically  scattered  light.  The  actual  time-variation  o 
droplet  diameter  will  be  a  function  of  the  droplet  nominal  diameter,  the  droplet  mdex  of 
refraction  (including  absorption),  the  modulation  frequency,  the  laser  intensity,  and  the  droplet 
liquid  and  surrounding  medium  thermal  and  mass  transfer  properties.  The  work  reported  here 
provides  theoretical  modeling  of  the  PMMS  procedure  that  wUl  a.)  assist  m  developmg  the 
relationship  between  the  measured  modulated  elastically  scattered  light  and  the  aerosol  cloud 
properties,  b.)  assist  in  optimizing  PMMS  design  parameters,  and  c.)  assist  in  verifymg  physica 
understanding  through  the  comparison  of  theoretical  predictions  with  correspondmg  controlled 
experimental  measurements. 

Section  II.  describes  the  theoretical  thermal  and  mass  transfer  modeling  for  a  single 
droplet  heated  by  an  intensity  modulated  laser  beam.  In  Sec.  III.,  the  elastic  hght  scattering  from 
a  polydispersed  aerosol  cloud  with  modulated  heating  is  considered.  A  set  of  demonstration 
calculations  are  presented  and  examined  in  Sec.  IV.  and  a  discussion  of  possible  foture  work  is 
given  in  Sec.  V.  Section  VI.  contains  a  list  of  references. 


2 


II.  Single  Particle  Modeling 

The  following  basic  assumptions  are  made. 

1 . )  The  droplet  is  spherical  and  consists  of  a  solution  of  a  nonvolatile  solute  (e.g.,  NaCl)  in  a 

volatile  solvent  (e.g.,  water).  a-* 

2. )  The  surrounding  medium  is  moist  air  with  time-invariant  temperature  ^d  relative  humidity. 

3. )  Plane  wave  illumination  on  the  droplet  by  the  infrared  (heating)  and  visible  (scattering)  laser 

beams  is  assumed. 

4. )  The  droplet  has  a  single  uniform  temperature.  . 

5. )  The  heat  and  mass  transfer  processes  at  the  surface  of  the  droplet  are  descnbed  usmg  a  quasi- 

steady,  boundary  layer  model.  ,  ,  ,  f 

6. )  A  form  of  Raoult's  law  is  used  to  determine  the  vapor  pressure  at  the  droplet  s  surface. 

The  notation  is  as  follows. 

Cd  =  specific  heat  per  mass  of  droplet  solution 
D  =  2  Tdo ,  droplet  initial  diameter 

=  log-normal  distribution  mean  diameter 
D  =  mass  diffusion  coefficient  for  vapor  in  air 
E  =  electric  field  vector  of  elastic  scattered  light 
hfg  =  enthalpy  of  evaporation  of  solvent 
H  =  magnetic  field  vector  of  elastic  scattered  light 
labs  =  time  dependent  intensity  of  the  infrared  (heating)  laser  beam 
Isca  =  intensity  of  the  visible  (elastic  scattering)  laser  beam 
k  =  thermal  conductivity  of  moist  air 
md=  (mw  +  ms)  =  droplet  mass 
nis  =  niso  ~  mass  of  solute  in  droplet  (time-invariant) 
mw  =  mass  of  solvent  in  droplet 
mwo  =  initial  mass  of  solvent  in  droplet 
Mw  =  molecular  weight  of  solvent 
Ms  =  molecular  weight  of  solute 

m"  =  mass  flux  of  vapor  at  the  droplet  surface  (positive  outward) 

miiq  =  initial  mass  of  aerosol  liquid  in  elastic  scattering  probe  volume 

N  =  number  of  droplets  in  elastic  scattering  probe  volume 

nabs  =  complex  index  of  refraction  of  droplet  at  the  infrared  wavelength 

n,ca  =  complex  index  of  refraction  of  droplet  at  the  visible  wavelength 

p(D)  =  probability  density  function  of  initial  aerosol  size  distribution 

Pg  =  saturation  pressure  of  solvent 

pv,d  =  vapor  pressure  at  droplet  surface 

Pv  00  =  vapor  pressure  at  ambient  condition 

Poo  =  ambient  pressure 

Qabs  =  droplet  absorption  efficiency  at  infrared  (heating)  laser  wavelength 
=  heat  flux  due  to  conduction  at  the  droplet  surface  (positive  outward) 
r<i  =  droplet  radius 
tdo  =  initial  droplet  radius 
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Rv  =  constant  for  vapor 
Td  =  droplet  temperature 
Tdo  =  initial  droplet  temperature 

Tf  =  film  temperature  of  the  boundary  layer  above  the  droplet  surface  =  (Td  +  Too )  /  2 
Too  =  temperature  of  the  ambient  air 
Vjo  =  average  droplet  volume  of  initial  aerosol  distribution 
Kbs  =  wavelength  of  infrared  (heating)  laser  beam 
Xsca  =  wavelength  of  visible  (scattering)  laser  beam 
^dir  =  polar  angle  direction  of  observed  elastic  scattered  light 
=  mass  density  of  droplet  solution 

mass  density  of  vapor  at  the  droplet  surface 
Pvoo~  density  of  vapor  in  the  ambient  air 

aoL  =  log-normal  distribution  standard  deviation  parameter 
<f>oo  =  relative  humidity  of  the  vapor  in  the  ambient  air 
<t)dir  =  azimuthal  angle  direction  of  observed  elastic  scattered  light 


An  application  of  conservation  of  energy  for  the  droplet  provides  the  following  equation 
for  the  time  rate  of  change  of  the  droplet  temperature, 

^=—^[hbs{^rl)Q,i,,  -  rh:i47rrl)hf,  -  ]  (1) 

dt  {rrid  Cd)  ^ 

where  the  three  terms  on  the  right-hand-side  represent,  respectively,  the  rate  of  laser  heating,  the 
rate  of  evaporative  cooling,  and  the  rate  of  conductive  cooling.  In  addition,  the  time  rate  of 
change  of  the  mass  of  the  solvent  of  droplet  is  the  negative  of  the  rate  of  solvent  evaporation. 


The  heat  flux  and  mass  flux  at  the  droplet  surface  are  evaluated  using  a  quasi-steady  boundary 
layer  assumption  so  that 


and 


( 


—  )^{Pv;d  ~  Pv,oo) 
Td 


(3) 


where  in  Eq.  (3), 

Pv,oo 

P'><°°  ~  ry  rp 
oo 

and 


'(5) 
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Pv,d 
~  R^Td 

(6) 

where 

Pv^oc  ”  4^ooPgi^oc) 

(7) 

and 

Pv,d  “  [1  ^ 

(8) 

Equation  (8)  is  a  representation  of  Raoult's  law,  where  the  assumption  is  made  that  the  vapor 
pressure  of  the  solvent  at  the  droplet  surface  is  equal  to  the  saturation  pressure  of  the  vapor  at  the 
droplet  temperature  times  the  mole  fraction  of  the  solvent  within  the  droplet.  The  factor  of  2  in 
the  denominator  of  Eq.  (8)  occurs  with  the  assumption  that  the  solute  is  ionic,  and  thus  each 
molecule  of  solute  effectively  divides  into  2  separate  ion  molecules. 

Before  the  laser  heating  is  initiated,  it  is  assumed  that  the  droplet  is  in  equilibrium  (both 
thermally  and  with  regard  to  composition)  with  the  surroundings.  The  mitial  mass  ratio  of  solute 
to  solvent  within  the  droplet  is  thus  determined  by  equating  Eqs.  (7)  and  (8)  with  the  initial 
condition  that  Td~  Too, 


{rasolrriyio)  =  ^  —  {^s/ (9) 

2(p 

For  a  droplet  of  initial  radius,  rdo,  the  initial  mass  of  solvent  and  solute  then  can  be  found  in 
terms  of  this  ratio, 


(mso/m,co)] 


(10) 


and 


rrisa  =  mwo{mso /mvio)  (H) 

Finally,  since  the  total  mass  of  the  droplet  is  the  sum  of  the  mass  of  the  solute  (which  is  constant) 
and  the  mass  of  the  solvent  (which  varies  in  time),  if  the  mass  of  the  solvent  in  known,  then  the 
radius  of  the  droplet  can  be  determined  from 


rd 


{niw  +  TOao)] 

{An/3)pd  J 


(12) 


Beginning  with  a  droplet  of  known  radius,  Eqs.  (9)  -  (1 1)  are  used  to  find  the  initial  mass 
of  solvent  and  mass  of  solute  within  the  droplet.  Once  the  laser  heating  starts,  the  time  variations 
of  droplet  temperature,  size,  and  composition  are  determined  by  solving  Eqs.  (1)  and  (2)  [along 
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with  the  associated  Eqs.  (3)  -  (8)  and  (12)],  using  a  finite  difference  predictor-corrected  method 
to  "inarch"  along  in  time.  The  following  property  relationships  are  also  required: 

Crf  =  func  (Td),  hfg  =  func  (To),  D  =  func  {Tf,Poo)>  k  =  func  (Tf), 

=  func  (Td,  {ms/niw)),  Qabs  =  func  (r d,  Uabs))  =  func{\abs-i  (juslmw))  where 

Tf  =  (Td  +  Too) /2  is  the  boundary  layer  film  temperature. 
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ni.  Elastic  Scattering  from  a  Polydispersed  Aerosol 

The  nondimensionalized  far-field  elastic  light  scattering  intensity  in  the  {ddir^  ^dir  ) 
direction  (relative  to  the  propagation  direction  of  the  incident  elastic  scattering  laser  beam)  for  a 
single  droplet  can  be  determined  from  Lorenz-Mie  theory  and  is  given  by 

Sr  =  lim  {r/rdf{—)Re{E  X  it)r/Isca  =  func  {Odir,  4>dir,  U,  ^sca,  fisca)  (13) 

Stt 

rlrj—^oo 

where  n,ca  =  func  (A,ca,  (m,/m^)) ,  If  p(D)  is  the  initial  droplet  size  distribution  for  the 
aerosol,  then 


roc 

K.o=(t)/  (m) 

6  J-oo 

provides  the  average  droplet  volume  in  the  initial  distribution  and  the  number  of  droplets  in  the 
probe  volume  is  then  given  by 


N  = 


miig 

PdO^dO 


(15) 


where  is  the  initial  amount  of  aerosol  hquid  in  the  elastic  scattering  probe  volume.  An 
expression  for  the  total  elastic  scattering  for  the  aerosol  can  then  be  determined  by  integrating 
over  the  individual  droplet  scattering  as  follows 


=  rlSMD)dD 

Tfl/linj-sca 


miigisca  PdO^dO 

While  any  aerosol  size  distribution  may  be  assumed,  for  the  calculations  presented  in 
Sec.  IV.,  a  log-normal  size  distribution  was  used  where 


(16) 


p{D)  = 


1  f  [ln{D/D,,)]\ 

V  Kl  ^ 


(^DL 


■exp 


ln{Di/ Dref)  —  ln{D  /  Dref) 


and 


”01  =  {\ln{Dlbi)f} 


1/2 


(17) 


(18) 


(19) 
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where  the  overbar  indicates  a  size  distribution  averaged  quantity.  As  an  example.  Fig.  1  shows  a 
log-normal  distribution  for  an  aerosol  with  a  log-normal  mean  diameter  (  Di)  of  2.0  microns  and 
a  log-normal  distribution  standard  deviation  parameter  (  upz,  )  of  0.20. 


8 


IV.  Demonstration  Calculations 


A  FORTRAN  computer  program  was  written  with  an  algorithm  based  on  the  procedures 
outlined  in  Secs.  11.  and  IH.  A  set  of  demonstration  calculations  were  performed  assuming  salt 
water  droplets  (solvent  =>  water,  solute  =>  NaCl)  with  CO2  laser  heating  =  10.6  //m,  n^bs 
=  (1.179,  0.072)  ]  and  He-Ne  laser  elastic  scattering  [A^ca  =  0.633  //w,  n^ca  =  (1.332, 1.47  e-08)]. 
Pure  water  properties  were  used  for  Cd(T),  hfg(T),  k(T),  D(T,p),  and  Pg(T).  The  droplet  mass 
density,  pd[T,(ms/mw)],  was  evaluated  including  both  the  effects  of  temperature  and  salt 
concentration.  Curvefits  for  the  absorption  efficiency,  Qabs ,  and  the  nondimensionalized  far- 
field  scattering  intensity,  ,  as  a  function  of  particle  size  parameter,  a  =  2  7t  rd  /  A ,  were 

generated  using  a  previously  developed  computer  program  (Lorenz-Mie, plane  wave  on  a  sphere 
electromagnetic  interaction).  A  graphical  presentation  of  Qabs(oCabs)  and  5r(ctsca)  are  given, 
respectively,  in  Figs.  2  and  3. 

The  nominal  conditions  for  the  calculations  that  follow  are  Too  (ambient  temperature)  — 
20  °C,  Poo  (ambient  pressure)  =  101.3  kPa,  and  Tdo  (initial  droplet  temperature)  =  20  “C  .  A  log¬ 
normal  droplet  distribution  is  assumed  with  a  log-normal  standard  deviation  p^ameter  of  Gdl  - 
0.20  and  a  log-normal  mean  diameter,^^,  of  either  0.50,  1.00, 2.00,  or  4.00  microns.  The 
infrared  (heating)  laser  beam  is  chopped  (i.e.,  square  wave  with  an  intensity  varying  from  zero  to 
lo)  with  a  nominal  intensity  of  =  1.0  MW/m^ .  The  total  elastic  far-field  scattering  intensity 
expression,  Sr)tot ,  was  computed  (nominally)  for  a  geometrical  arrangement  corresponding  to 
pure  backscattering  (0dir  =180",  (j)dir  =  90" ). 

A.  Time  Constants 

Characteristic  time  constants  for  laser-on  ( Ton )  snd  for  laser-off  (  Toff )  were  determined 
by  using  a  low  frequency  modulation  (2.5  Hz  or  5.0  Hz)  so  that  the  aerosol  size  distribution 
would  reach  equilibrium  before  the  laser  is  switched  back  offron.  The  time  constant  was  defined 
as  the  time  it  took  the  far-field  scattering  to  recover  63.2  %  (i.e.,  1  -  e‘*  )  of  the  overall  difference 
of  the  laser-on  versus  laser-off  far-field  scattering.  Figure  4  shows  the  low  frequency  time 
variation  of  the  far-field  total  scattering  for  the  four  different  log-normal  mean  diameter  size 
distributions.  A  plot  of  the  corresponding  calculated  time  constants,  both  for  laser-on  and  for 
laser-off,  versus  mean  diameter,  is  given  in  Fig.  5.  As  seen  in  Figs.  4  and  5,  the  laser-on  (i.e., 
heating)  time  constant  increases  almost  linearly  with  increasing  aerosol  mean  diameter.  For  the 
gmall  diameters  (0.5  micron  and  1.0  micron),  the  laser-off  and  laser-on  time  constants  are 
approximately  equal.  However,  for  mean  diameters  greater  than  1.0  micron,  the  laser-off  time 
constant  is  significantly  longer  than  the  laser-on  time  constant  and  the  laser-off  time  constant 
may  either  increase  or  decrease  with  increasing  mean  diameter  (see  Fig.  5). 

The  time  constant  results  in  Figs.  4  and  5  were  for  a  nominal  infrared  (heating)  laser 
intensity  of  1 .0  MW/m^.  The  effect  of  laser  intensity  on  the  time  constants  was  investigated  ^d 
the  results  are  shown  in  Figs.  6  and  7.  For  the  1.0  micron  diameter  distribution  (Pig-  6),  the  tiine 
constants  decreased  slightly  with  increasing  laser  intensity.  For  the  2.0  micron  distribution  (Fig. 
7)  the  laser-on  time  constant  decreased  significantly  with  increasing  laser  intensity  while  the 
laser-off  time  constant  increased  slightly  with  increasing  laser  intensity. 
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The  time  constant  results  in  Figs.  4  and  5  were  for  a  nominal  ambient  relative  humidity  of 
90%.  The  effect  of  ambient  relative  humidity  on  the  time  constants  was  investigated  and  the 
results  are  shown  in  Figs.  8  and  9.  For  the  1.0  micron  diameter  distribution  (Fig.  8),  the  time 
constants  (both  laser-on  and  laser-off)  increased  significantly  with  increasing  ambient  relative 
humidity.  For  the  2.0  micron  distribution  (Fig.  9),  the  laser-on  tiine  constant  increased  slightly 
and  the  laser-off  time  constant  increased  significantly  with  increasing  ambient  relative  humidity. 


B.  Frequency  Response 

Systematic  calculations  were  also  performed  investigating  the  effects  of  laser  heating 
modulation  frequency  on  the  subsequent  total  far- field  scattering.  Figure  10  shows  the  0.5 
micron  mean  diameter  total  far-field  scattering  intensity  as  a  function  of  time  for  the  nominal 
conditions  and  modulation  frequencies  of  5, 10, 20, 50, 100,  and  200  Hz.  As  seen  in  Fig.  10.,  for 
the  0.5  micron  mean  diameter,  the  peak-to-peak  far-field  scattering  intensity  is  approximately 
constant  for  frequencies  up  to  50  Hz,  and  then  decreases  slightly  as  the  frequency  is  increased  to 
100  Hz  and  then  decreases  significantly  as  the  frequency  is  increased  to  200  Hz.  Simlar  time- 
varying  far-field  scattering  intensity  plots  are  provided  for  mean  diameters  of  1 .0  microns,  2.0 
micron,  and  4.0  micron,  respectively,  in  Figs.  11, 12,  and  13.  The  "cut-off"  frequency  at  which 
the  total  far-field  scattering  begins  to  decrease  apparently  decreases  with  increasmg  mean 
diameter.  This  observation  is  quantified  in  Fig.  14,  where  the  normalized  standard  deviation  of 
the  total  far-field  scattering  intensity  (normalized  relative  to  the  mean  value  of  the  total  far-field 
scattering  intensity)  is  plotted  versus  frequency.  As  can  also  be  seen  in  Figs.  12  and  14,  the 
relative  variation  of  total  fer-field  scattering  is  particularly  large  for  the  2.0  micron  mean 
diameter  distribution. 
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V.  Possible  Future  Work 


A  theoretical  modeling  of  the  photothermal  modulation  of  Mie  scattering  spectroscopy 
technique  has  been  developed  and  demonstration  calculations  have  been  performed.  Appropriate 
simplifying  assumptions  were  applied  in  developing  this  initial  modeling.  Possible  future  work 

includes  the  following.  j-  j  i 

1 . )  Direct  comparisons  between  ejq)erimental  measurements  and  correspondmg  model 
calculations  are  needed  to  verify  the  overall  correctness  of  the  model  and  to  verify  the  validity  of 
the  simplifying  assumptions. 

2. )  The  form  of  Raoult's  law  used  in  Sec.  II  is  valid  for  low  solute  concentrations.  A  thorough 
literature  search  could  be  performed  to  see  if  a  more  rigorous,  and  more  wide-ranging,  farm  of 
Raoult’s  law  can  be  identified.  In  addition,  the  effect  of  surface  tension  on  Raoult’s  law  was 
omitted  (probably  a  good  assumption  except  for  very  small  droplets)  and  this  effect  could  be 

included  in  a  more  advanced  form  of  Raoult’s  law.  •  j  ^ 

3. )  The  effects  of  solute  concentration  on  droplet  properties,  in  particular  the  complex  mdex  of 
refraction  at  both  the  laser  heating  and  laser  scattering  wavelengths,  needs  further  study. 

4. )  A  finite-difference  analysis  could  be  performed  to  investigate  the  correctness  of  the  quasi¬ 
steady  boundary  layer  assumption  for  the  modeling  of  the  mass  transfer  and  heat  transfer  at  the 
droplet  surface.  If  necessary,  corrections  could  be  incorporated  into  the  model. 

5. )  A  uniform  droplet  temperature  was  assumed.  This  is  probably  a  good  assumption  for 
relatively  small  droplets  and  low  modulation  frequencies.  However,  the  possible  effects  of 
distributed  heating  could  be  investigated. 
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Fig.  1.  Log-normal  probability  density  distribution  function  with  a  log-normal  mean  diameter 
(Di)  of  2.0  microns  and  a  log-normal  standard  deviation  parameter  (  o^dl  )  of  0.20. 
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Fig.  2.  Absorption  efficiency  as  a  function  of  size  parameter,  X^bs  =  10.6  pm, 
(1.179,0.072). 


Fig.  3.  Nondimensionalized  far-field  scattering  intensity  (S^)  for  pure  backscatter  (Ga 
180°,(|)dir  =  90°)  as  a  function  of  size  parameter,  Xsca  =  0.633  pm, 

=  (1.332,1.47E-08). 


Fig.  6.  Laser-on  and  laser-off  far-field  scattering  time  constants  as  a  function  of  heating  laser 
intensity.  Log-normal  mean  diameter  equal  to  1 .0  micron. 


Fig.  7.  Laser-on  and  laser-ofif  far-field  scattering  time  constants  as  a  function  of  heating  laser 
intensity.  Log-normal  mean  diameter  equal  to  2.0  micron. 


Fig.  9.  Laser-on  and  laser-off  far-field  scattering  time  constants  as  a  function  of  ambient  relative 
humidity.  Log-normal  mean  diameter  equal  to  2.0  micron. 
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Fig.  10.  Total  far-field  scattering  versus  time  for  modulated  heating  laser  intensity  at  modulation 
frequencies  of  5,  10, 20,  50,  100,  and  200  Hz.  Log-normal  mean  diameter  equal  to  0.5  micron. 


Fig.  12.  Total  far-field  scattering  versus  time  for  modulated  heating  laser  intensity  at  modulation 
frequencies  of  5, 10, 20,  50, 100,  and  200  Hz.  Log-normal  mean  diameter  equal  to  2.0  micron. 
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Fig.  13.  Total  far-field  scattering  versus  time  for  modulated  heating  laser  intensity  at  modulation 
frequencies  of  5, 10, 20,  50, 100,  and  200  Hz.  Log-normal  mean  diameter  equal  to  4.0  micron. 
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Fig.  14.  Normalized  standard  deviation  of  total  far-field  scattering  intensity  as  a  function  of 
modulation  frequency  for  log-normal  mean  diameters  of  0.5, 1.0, 2.0,  and  4.0  microns. 


